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INTRODUCTION 


Hxcited  states  of  interacting  systems,  that  are  bound  in  the  ground  state 
only  by  van  der  Waals  forces,  can  interact  to  form  strongly  bound  species.  This 
is  found  even  for  systems  in  which  only  one  of  the  fragments  is  excited,  and  is 
exemplified  by  the  exciiner  states  He-He*,  Ar-Ar*.  and  IT, -I  F*.  Here  the  asterisk 
denotes  the  first  excited  state  of  the  same  spin  symmetry  as  the  ground  state. 
Interest  in  the  latter  system  has  increased  recently  with  the  study  of  Nicolaides, 
Theodorakopoulos,  and  IYtsalakis  (NTP)1  of  the  1F(X  'XR')  1F(B  *iJu4)  system 
because  of  the  ionic  character  of  the  B  state  which  these  authors  label  at  -1.0a0. 
where  charge  transfer  occurs,  a  maximum  ionieity  excited  states  (MHOS);  see  Fig. 
1  for  the  correlation  diagram.  This  state  is  of  special  interest  because  of  its 
strong  electrostatic  binding. 

A  model1  3  based  on  MHOS  properties  suggests  that  bound  excited  states  of 
polyatomic  systems  can  be  formed  in  regions  characterized  by  an  avoided  crossing 
with  the  ground  state,  if  one  of  the  interacting  molecules  can  exist  in  a  MHOS. 
The  description  for  IF  IF*,  where  IF*  denotes  IF(B  ’X,/  )  is  one  in  which  a  posi¬ 
tive  ion  complex  II3!  is  formed,  and  interacts  in  its  ground  state  equilibrium 
geometry  (equilateral  triangle  with  r  ~  1.6oa0)  with  II  .  The  MIFS  geometry 
corresponding  to  Il+  H  for  the  parent  H2(B  ‘E,/  )  is  the  charge  transfer  region  at 
4.()a0.  The  mechanism  supported  by  NTP  with  Cl  computations  is  one  of  I lj+ , 
which  is  electron  deficient  in  the  center  of  the  triangle,  interacting  with  H  at  a 
distance  of  roughly  1.0a,,  above  the  plane  of  the  ||{’  triangle. 

] 


During  the  contract  period  studies  of  the  I L*  I h  system  have  been  pursued 
in  this  laboratory  using  the  lixed-node  quantum  Monte  Carlo  (FNQMC)  method. 
A  brief  introduction  to  the  method  is  presented  and  followed  by  a  description  of 
FNQMC  and  other  calculations  of  the  subject  system. 

QUANTUM  MONTE  CARLO 

Monte  Carlo  approaches  to  solving  problems  with  many  degrees  of  freedom 
are  a  class  of  statistical  methods  having  in  common  the  generation  of  “random" 
numbers.  In  the  past  few  years,  Monte  Carlo  approaches  have*  seen  increased 
application  in  a  number  of  diverse  fields.  What  we  mean  here  by  quantum 
Monte  Carlo  (QMC)  is  a  Motile  Carlo  procedure  which  solves  the  Sehrodinger 
equation.  'Phis  is  to  be  distinguished  from  so-called  variational  Monte  Carlo,  in 
which  one  obtains  expectation  values  for  a  given  trial  wave  function. 

This  ability  to  stochastically  solve  the  Sehrodinger  equation  provides  an 
alternative  to  conventional  techniques  of  quantum  chemistry.  Karly  work1  has 
shown  that  highly  accurate  total  energies  and  correlation  energies  can  be 
obtained  by  QMC.  In  fact,  in  a  proeedu rally  simple  manner,  accuracies  exceeding 
those  of  the  best  ah  initio  configuration  interaction  calculations  have  been 
obtained. 

The  essence  of  the  procedure  is  to  simulate  a  quantum  system  by  allowing  it 
(and  an  ensemble  of  differently  prepared  systems)  to  evolve  under  the  time- 
dependent  Sehrodinger  equation  in  imaginary  time.  Lxcited  slat 


es  of  I  lie  same 


By  writing  tin-  imaginary-time  Schrodinger  equation  with  a  shift  in  the  zero 


of  energy  as 

]  =  DvXliJ  )  +  [lvr  \'(!i)]^([i.T).  ( I  ) 

we  see  that  it  may  lie  interpreted  as  a  generalized  diffusion  equation.  The  lirst 
term  on  the  right-hand  side  is  the  ordinary  diffusion  term,  while  the  second  term 
is  a  position-dependent  rale  (or  branching)  term.  Bor  an  electronic  system. 
I)  --  fr’/Sni,,  !i  is  the  three-N  dimensional  coordinate  vector  of  the  N  electrons, 
and  \  (H)  is  tin'  Coiilomh  potential.  Since  diffusion  is  the  continuum  limit  of  a 
random  walk,  one  may  simulate  Fvp  (l)  with  the  function  T  (note,  not  T")  as  the 
density  oi  "walks  .  The  walks  undergo  an  exponential  birth  and  death  as  given 
by  the  rate  term. 

I  he  steady-state  solution  to  Ivp  (I)  is  the  time-independent  Selirddinger 
ecpiation.  Thus  we  have  'l'(fLt)— *■  <^>( LL)»  where  (t>  is  an  energy  eigenstate.  The 
value  of  Ivp  at  which  the  population  of  walkers  is  asymptotically  constant  gives 
the  energy  eigenvalue.  Karly  calculations  employing  Ivp  (I)  in  this  way  wen 
done  by  Anderson  on  a  number  of  one-  to  four-elect  roll  systems.' 

lu  order  to  treat  systems  larger  than  two  eleeinms,  tin  fermi  nature  (it  the 
electrons  must  be  taken  into  account.  Ihe  antisymimt  r\  of  the  eigenfunction 
implies  that  T  must  change  sign:  however,  a  density  (e.g..  of  walkers)  cannot  be 
negative.  Ihe  method  which  imposes  antisymmetry,  and  at  the  same  time  pro¬ 
vides  efficient  sampling  (thereby  reducing  llm  statistical  "noise").  B  importance 
sampling  with  an  antisymmetric  trial  function  T|.  The  zeroes  (nodes)  of  \f». 
become  absorbing  boundaries  for  1  he  d i 1 1 1 1  ion  piece  .  which  maintains  i  l,e 


«.  .1. 


MM  IJ’*.!'! 


^.1  J.F  A.M.MJ  t.1 


antisvmmetrv.  Tin*  :u  1  <  11 1  ion  a  I  boundary  condil  ion  that  'P  vanish  :it  the  nodes  of 


'I'-p  is  I  In*  lixed-node  :  i  p  p  r<  >  x  i  1 1 1 :  >  t  ion.  The  magnitude  of  1 1 1  <  •  error  thus  ini  rod  need 


depends  on  I  lie  accuracy  of  the  nodes  of  4».p(R),  and  vanishes  as  T  |  approaches 


the  time  eigenfunction.  To  the  extent  that  'f'-p  is  a  good  approximation  of  the 


wave  function.  I  lie  true  eigenfunction  is  almost  certainly  'piite  small  near  the 


nodes  of  4>t.  Thus  one  expects  the  lixed-node  error  to  l>e  small  for  reasonable 


choices  of  4'  p. 


To  implement  importance  sampling  and  the  lixed-node  approximation.  Imp 


(I)  is  multiplied  on  both  sides  by  T-p,  and  rewritten  in  terms  of  t  he  new  probabil¬ 


ity  density  f(R,t)  =  vl'-p(R )T(R,t).  The  resultant  equation  for  f(R.t)  may  be  writ- 


|L  =  I)v-f+[IOT-  h:L(R)]f-  Dv|fl-’q(K)i. 


Phe  local  energy  Rp(R),  and  the  “quantum  force’ 


’q(R)  tire  simple  functions  of 


'I'-p  given  by 


I'lXR)  =  lN'T(!i)/*T(ii). 


LI)  =  2v*t(R-)/'1't(R)-  (3b) 

hqualion  (2).  like  Imp  ( I  )  is  :i  generalized  <  li  1 1  u  si  <  >  1 1  equation,  though  now  with  I  lie 


addition  of  a  drift  term  due  to  the  presence  of  f'0. 


order  to  perform  the  random  walk  implied  by  Imp  (2)  we  use  a  short-time 


approximation  to  the  (ireeu’s  function  which  is  use. I  to  evolve  f(R,t)  - 1'( R /  .1  (  ,*). 


I  his  •■volution  process  is  iterated  to  large  |.  The  (ireeu  s  function  lieeomes  exact 


in  the  limit  ol  vanishim:  timc-->tep  si/e. 


X 


1.075 


for  C  symmetry  are  indicated  in  parentheses.  The  lack  of  coincidence 
the  degenerate  E  curves  reflects  their  calculation  in  lower  symmetry  a 
provides  an  indication  of  Mf.SC.F  oovergence. 


Geometries  and  energies  for  the  approach  of  Ho(B)  to  H2(X) 


step 

R/ 

Li  L, 

Lo(+) 

m-) 

D 

energy  lowing 

I 

5-2.2 

1.40  2.43 

1.215 

1.215 

0.0 

-0.01 

n 

2.2 

1.40  2.43 

1.215-1.03 

1.215-0.5 

0.715 

-0.20 

m 

2.2 

1.60  2.63 

1.03-2.13 

0.5 

0.865 

-0.20 

IV 

2.2-1.75 

1.80  2.63 

2.14 

0.5 

0.865 

-0.20 

V 

1.75-1.55 

1.60  2.73 

2.23 

0.5 

0.015 

-0.48 

VI 

1.55 

1.60  3.13 

2.23-2.63 

0.5 

1.115 

-0.69 

vn 

1.50 

1.70  3.23 

2.63-2.73 

0.5 

1.165 

-0.11 

TOTAL  ENERGY  LOWERING  FOR  THE  SEVEN  STEPS  IS  5.75 
ev. 


•  Distances  in  a.u.;  energies  in  ev. 

R'  -  distance  between  the  midpoint  of  H2(X)  and  the  point  where  Hj(B)  meets 


the  H2(X)  plane. 

Lt  -  length  of  Hj(X). 

L,  •  length  of  H^B). 

L,(+)  -  length  of  Hj(B)  above  Hj(X)  plane. 

L^f-)  -  length  of  H}(B)  below  H^(X)  plane. 

D  -  magnitude  of  shift  of  HjfB)  midpoint  (  above  (+)/below(-)  )  plane  of  HjQC). 


asymptotic  region  to  the  neighborhood  of  the  H4  MIES  configuration.  The  table 
also  contains  the  energy  lowering  associated  with  the  steps  that  are  plotted  in 
Fig.  6.  It  is  noteworthy  that  no  energy  barrier  is  encountered  along  this  path. 
Further,  Fig.  7  show  that  when  H(B)  is  displaced  from  bisecting  H2(X)  in  step  II, 
charge  transfer  immediately  occurs.  (Similar  behavior  has  been  found  for  the 
related  H2(B)  +  He  system  in  an  independent  study  by  WAL.) 

FNQMC  STUDY  OF  THE  GROUND-STATE  PYRAMIDAL 
STRUCTURE  FOR  SYMMETRY 

FNQMC  calculations  using  the  MCSCF  trial  functions  discussed  above  yield 
0.6-1.0  eV  energy  lowering  compared  to  the  results  of  NTP  and  are  presented 
in  Fig.  8.  Such  a  large  change  was  not  anticipated  and  so  it  was  important  to 
test  the  validity  of  this  finding.  To  this  end  a  configuration  interaction  calcula¬ 
tion  including  all  single  and  double  excitations  (SDCI)  using  the  MCSCF  pilot 
study  basis  set  was  carried  out  at  R  =  3.4  a.u.  The  energy  was  0.32  eV  lower 
than  NTP’s  value  and  is  consistent  with  the  improvement  expected  based  on  stu¬ 
dies  of  other  systems.  The  FNQMC  results  of  Fig.  8  obtained  using  a  new  trial 
function  optimization  algorithm,  mentioned  in  the  next  section,  are  generally  an 
improvement  over  those  of  Fig.  8  obtained  using  MCSCF  trial  functions. 

FNQMC  STUDY  OF  THE  EXCITED  STATE 


These  calculations  provide  the  severest  test  of  the  FNQMC  approach 
because  of  the  lack  of  knowledge  of  the  accuracy  of  the  excited  state  trial  func¬ 
tion  needed  to  provide  a  nodal  description  that  assures  orthogonality  to  the 
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ground  stale  of  tin-  same  (A1  )  symmetry  (C  ).  MCSCF  calculations  close  to  tin- 
avoided  crossing  sullered  from  root  flippl  11^4.  De-pitc  (lie  use  of  familiar  MCSCF 
strategies  to  address  the  problem,  it  could  not  be  resolved. 

The  MCSCF  convergence  problem  had  been  encountered  earlier  for 
He  +  ||,(U)  but  resolved  by  the  use  of  an  ab  initio  Cl  method.  Interest  here  in 
using  the  FNQMC  method  led  us  to  develop  a  method  for  trial  function  construc¬ 
tion6  that  avoids  the  MCSCT  procedure.  The  approach  amounts  to  the  introduc¬ 
tion  of  parameter  optimization  in  the  random  walk  process.  Using  group  theory, 
a  projection  operator  is  constructed  and  used  to  constrain  the  wave  function  to 
have  the  symmetry  properties  of  the  state  of  interest.  This  method  has  been 
employed  in  the  present  effort  to  generate  the  A-state  FNQMC  results,  denoted 
A(FNQMC),  of  Fig.  8.  Calculations  using  this  method  are  in  progress  to  com¬ 
plete  the  E-state  curve,  labeled  E(FNQMC),  in  Fig.  8  for  3.S  <  R  <  6.0  a.u. 

SA-MCSCF  STUDY  OF  DISTORTED  GEOMETRIES 

Calculations  have  been  carried  out  using  the  state  averaged(SA)-MCSCF 
method  to  develop  trial  functions  for  a  QMC  study  of  the  topography  of  the 
ground-  and  excited-state  potential  energy  surface  (pes)  in  the  region  of  the  sym¬ 
metric  geometry  of  the  MIES  determined  by  NTP.  Figure  <J  presents  the  coordi¬ 
nate  system.  Fig.  10  provides  perspective  views,  and  Fig.  1 1  (R  =  4.0  a.u.)  and 
Fig.  12  (H  —  3.8  a.u.)  show  contour  maps  of  the  pcs  in  the  MIES  region.  Fig¬ 
ures  11  and  12  show  that  the  ground  state  has  a  saddle  point  at  smaller  R  than 
the  minimum  of  the  excited  state  and  that  both  features  correspond  to  an  isos¬ 
celes  triangle  base  for  the  MIES  system.  Further  geometry  optimization  is 


3-D  POTENTIAL  ENERGY  SURFACE  (R=3.8  a.u.) 


wr\ri».’irir 


P-  u-t  I 


I  7 


Afc^fc 


**  A 


UU 


CONTOUR  MAP  (R— 4.0  a.u.) 


r- — t-t-— 

W  w  41  A  1/  W  *k  4*4-rU  ^iV“ 

1  ! 

/s/jfflm 


's' 


7  -» 

if* 


— ^  \ 

■«•  \  >  / 

\\ 

•  >  i  // 
/  i  //> 


:3I1I! 


''&Mw 

.'//AMM 


w/Mfi 

/  /////,®M® 


explored  in  I  ig.  Ft  which  plots  p<  tented  energy  as  a  In  mt  i<  >n  of  <1.  tin  •  li  -  j  ■ ! :  <  ■  ■< 
llifllt  ( -f  li,  ill*-  I'll.-'  of  (In-  i|(  is,,-.-.  !t  .  arra  nge  m »  lit.  Tin-  1 1  ■  i  1 1  i  1 1 1 1 1 1 :  i  f -r 

the  excited  stale  is  found  for  .1  —  0.1  :i.ii. 

NONADIABATIC  COUPLING 

'I  lie  stability  of  the  MILS  system  is  dependeii!  <  n  iionadiabatie  couplim' 
(NA(')  to  the  ground  slate.  Although  the  locus  of  this  study  is  ehar.'e't  eri/at  ],  a 
of  the  re-o ion  of  pes  of  the  MILS,  the  need  to  ascertain  lil'st  the  stability  of  tie 
excited  system  is  a  1 1  i  .e,  1 1  priority  here  hecau— ■  of  the  eflort  asso<-iat  m|  with  t  h< 
trial  function  optimi/.ation  method  and  tin  expense  of  OM( '  calculations  with 
small  statistical  variances.  To  calculate  tin-  NA<  '  matrix  element-  w  tool-, 
advaiit ace  of  the  simplicatiotis  made  possible  by  the  use  of  normal  mode  distor¬ 
tions  (Lie.  |.|)  from  the  highly  symmetric  (' M  symmetry.  Our  calculation  show 
that  nuclear  displacements  associatid  with  modes  Q,  and  Q-.  lead  to  strong  cou¬ 
pling  between  the  A  and  L  states.  Relatively  weaker  coupling  is  found  in  Q ,  and 
Qv  The  two  totally  symmetric  modes  ((.),  and  give  ie>  contribution. 

The  size  of  I  he  .\AC  matrix  elements  connecting  the  A'  state  to  the  L-tat. 
components  (Cu  notation  is  used  lu  re  to  indicate  parent  ,'c|  are  tahnlaled  in 
1  io, .  la  for  the  minimum  energy  geometry  of  the  A  state.  Tlu-..  results  e-talli  h 
that  there  is  strong  coupling  hetween  these  states.  Further  computational  study 
is  needed  to  confirm  these  predictions  and  to  estimate  reliably  the  excited  state 


lift  t  ime. 


<e2  A'  > 


-2.01 12(64) 
0.0017(0) 
0.0000 
-0.3901(389) 
-0.0039(13) 
0.0000 


A'  > 


0.0015(0) 
2.0127(98) 
0.0000 
-0.0010(3) 
-0.31 13(196) 
0.0000 


10 1  and  E2  arc  doubly  degenerate  states  where  E(  is  sym 
metric  and  U2  is  antisymmetric. 


Fi^uro  l!>.  ('1  font r i but  ion  to  coupling  mntrix  elements  For  normal  modi 

d i spl ncement  4  see  Fig.  14. 
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